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Abstract 

The Dirac quasinormal modes (QNMs) of the Kerr-Newman black hole are investigated using 
continued fraction approach. It is shown that the quasinormal frequencies in the complex uj plane 
move counterclockwise as the charge or angular momentum per unit mass of the black hole increases. 
They get a spiral-like shape, moving out of their Schwarzschild or Reissner-Nordstrom values and 
"looping in" towards some limiting frequencies as the charge and angular momentum per unit mass 
tend to their extremal values. The number of the spirals increases as the overtone number increases 
but decreases as the angular quantum number increases. It is also found that both the real and 
imaginary parts are oscillatory functions of the angular momentum per unit mass, and the oscillation 
becomes faster as the overtone number increases but slower as the angular quantum number increases. 
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I. INTRODUCTION 



It is well known that QNMs possess discrete spectra of complex characteristic frequencies 
which are entirely fixed by the structure of the background spacetime and irrelevant of the 
initial perturbations Q] 3|. Thus, it is believed that one can directly identify the existence of a 
black hole by comparing QNMs with the gravitational waves observed in the universe, as well 
as test the stability of the event horizon against small perturbations. Meanwhile, it is generally 
believed that the study of QNMs may lead to a deeper understanding of the thermodynamic 
properties of black holes in loop quantum gravity ^| j^f since the real part of quasinormal 
frequencies with a large imaginary part for the scalar field in the Schwarzschild black hole is 
equal to the Barbero-Immirzi parameter , a factor introduced by hand in order that 

the loop quantum gravity reproduces correct entropy of the black hole. It was also argued that 
the QNMs of Anti-de Sitter (AdS) black holes have a direct interpretation in terms of the dual 
conformal field theory (CFT) 0, 0, 0|. According to the AdS/CFT correspondence, a large 
static black hole in asymptotically AdS spacetime corresponds to a (approximately) thermal 
state in CFT, and the decay of the test field in the black-hole spacetime corresponds to the 
decay of the perturbed state in CFT. The dynamical timescale for the return to the thermal 
equilibrium can be done in the AdS spacetime, and then translated into the CFT using the 
AdS/CFT correspondence. Thus, QNMs for various black holes have been studied extensively 

Q-Q. 

Although many authors studied the QNMs of scalar, electromagnetic and gravitational per- 
turbations, the investigation of the Dirac QNMs is very limited [30, 31. Further- 
more, in these papers the study of the Dirac QNMs was limited to use Poshl- Teller potential 
or WKB approximative methods. The Poshl- Teller potential approximative may be a valuable 
tool in certain situations, but in general results obtained from the approximating potentials 
should not be expected to be accurate. The WKB method can be generalized to deal with 
the Reissner-Nordstrom and Kerr black holes. In the latter case, unfortunately, the method is 
not very simple to apply and does not yield very accurate results. It is well known that the 
Leaver's continued fraction technique is the best workhorse method to compute highly damped 
modes and is very reliable. However, the method can not be directly used to study the QNMs 
of the Dirac fields if we take the standard wave equations. The reason to obstruct the study of 
the Dirac QNMs using continued fraction technique is that for static black holes the standard 
wave equation 

d 2 



possesses a special potential 



(h3 +u*)Z± = V±Z ± (1.1; 



V ± = \ 2 -±\—\- 1.2 

r 4 cir* r z 
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which is the function of where A is a function related to the metric of the 

spacetimes, say A = r 2 — 2Mr for the Schwarzschild black hole. We can not investigate the 
Dirac QNMs using the continued fraction approaches because we have to expand the potential 
as a series at the event horizon but the factor vA does not permit us to do that. 



Recently, we found that the wave function and the potential of the Dirac fields can 







be expressed as new forms. Starting from the new wave function and potential, we J 
studied the Dirac QNMs of the Schwarzschild black hole using the continued fraction |38|] and 
HiU-determinant approaches 00 aad the Dirac QNMs of the Schwarzschild-anti-de Sitter 
and Reissner-Nordstrom-anti-de Sitter black holes using Horowitz-Hubeny approach 11]. The 
main purpose of this paper is to extend the investigation to the Kerr-Newman black hole 
because this black hole is the only asymptotically flat electro-vacuum solution of the Einstein- 
Maxwell system 4l| and it is the most general black hole among classical black holes. 

The organization of this paper is as follows. In Sec. 2 the decoupled Dirac equations and 
corresponding wave equations in the Kerr-Newman spacetime are obtained using Newman- 
Penrose formalism. In Sec. 3 the angular and radial continued fraction equations are intro- 
duced. In Sec. 4 the numerical results for the Dirac QNMs in the Reissner-Nordstrom, Kerr 
and Kerr-Newman black holes are presented. The last section is devoted to a summary. 



II. DIRAC EQUATIONS IN THE KERR-NEWMAN SPACETIME 

In a curve spacetime the Dirac equations [42] can be expressed as 

V2V B b'P B + ifjQB' = 0, 

V2V BB ,Q B + i/iP B > = 0, (2.1) 

where Vbb' is covariant differentiation, P B and Q B are the two-component spinors representing 
the wave functions, P B > is the complex conjugate of Pb, and \i is the mass of the Dirac particle. 
In the Newman-Penrose formalism [43] the equations become 



(D + e- 


-p)P° + (5 + it- 


cOP 1 = 




(A + /i 


- 7 )P 1 + (5 + /3- 


-r)P° 


= ^-^vq ' 


(D + e- 


- p)Q°' + (6 + 7T - 


- a)Q 1 ' 


= -2- 1 / 2 z^P 1 


(A + /Z 


-l)Q V + + P 


- r)Q ' 


= 2- l ' 2 ipP Q . 
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The null tetrad for the Kerr-Newman black hole can be taken as 

i»=( r2 + a2 1 o -] 

V A ' ' U ' A J ' 

n^ = -^((r 2 + a 2 ), -A, 0, a), 
2pp 



m * = ^( iasM , , (2.3) 



with 



p = r — iacos9, 

A = r 2 - 2Mr + a 2 + Q 2 = (r - r + )(r - r_), (2.4) 

where M, Q and a represent the mass, charge and angular momentum per unit mass of the 
Kerr-Newman black hole respectively, and r± = M ± \J M 2 — a 2 — Q 2 . As Q — > the Kerr- 
Newman metric reduces to the rotating Kerr metric, and as a — > it reduces to the charged 
Reissner-Nordstrom metric. Then the wave functions can be taken as 



P° = -M_ 1/2 (r)S_ 1/2 (^e- i( ^- m¥3) , 

P 

P 1 = R +1/2 (r)S +1/2 (6)e~^ t - m *\ 
Q 1 ' =R +1/2 (r)S_ 1/2 (6)e-^ t - m *\ 

go' = -iR_ 1/2 ( r )5 +1/2 (0)e- l ^- m ^, (2.5) 
p 

where u and m are the energy and angular momentum of the Dirac particle. After the tedious 
calculation Eq. (|2.2j) can be simplified as 

y/AV R- 1/2 = (A + z>r)VAM+i /2 , (2.6) 

VApJ(VAM +1/2 ) = (A - z/ir)R_ 1/2 , (2.7) 

A/2S+1/2 = -(A - aiMcos6)S-i/ 2 , (2.8) 

£{/ 2 S'_i/2 = (A + aficos9)S +1 / 2 . (2.9) 



with 



9 ziC n dA 

n <9r A A dr ' 

x d iK n dA 

Tr = 1 1 

n dr A A dr ' 

£ n = — — aojsinO H + n cot 

ad sin 6* 

$n = "^t + ausinO : — - + n cot 0, 

ao 1 sin v 

K = (r 2 + a 2 )uo — ma. (2.10) 
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We can eliminate S+1/2 (or S-1/2) from Eqs. ()2.8|) and (|2.9|) . Defining w = cos#, we find 
that the angular equation can be expressed as 

\2" 



d 

du 



dS s 
du 



+ 



(cluju) — 2ausu + s + A 



Irn 



[m + su) 
1 — u 2 



S, = 0, 



(2.11] 



where Ai m is the angular separation constant. It is interesting to note that the angular equation 
(12.11}) is the same as in the Kerr case. 

We will focus our attention on the massless Dirac field in this paper. Therefore, we can 



eliminate 



-1/2 



or vAM + i/ 2 ) from Eqs. (|2.fi|) and ()2.7j) to obtain a radial decoupled Dirac 



equation for vAl + i/2 (or R_i/2). Then, introducing an usual tortoise coordinate 



dr* 



and resolving the equation in the form 



r 2 + a 2 
A 



dr. 



(2.12) 



A -s/2 

\/r 2 + a 2 



(2.13) 



we obtain the wave equation 



d 2 V s 
dr 2 



+ (u 2 - Vs)V 8 = 0, 



with 



dZ . A 

L £ z _| 

dr* (r 2 + a 2 ) 2 



2mau(r 2 + a 2 ) — m 2 a 2 isK dA 
A h A dr 



Aisur + 



(2.14) 



(2.15) 



where Z = ~ 2 (r 2 +a 2 ) lr — (r 2 +a 2 ) 2 an< ^ ^ 2 = ~ ^ mauJ + (cl^) 2 ■ We will see that we can 
easily work out the Dirac quasinormal frequencies of the Kerr-Newman black hole from Eqs. 
(I2.11j) and ()2.14|) using the continued fraction approach. 



III. ANGULAR AND RADIAL CONTINUED FRACTION EQUATIONS 

For the Kerr-Newman black hole, we have to solve Eq. ()2.11|) numerically following Leaver 
Boundary conditions for Eq. ()2.11|) are that S s is regular at the regular singular points 
u = ±1 where the indices are given by ±(m + s)/2 at u = 1 and ±(m — s)/2 at u = —1. So a 
solution to Eq. (J2.11j) may be expressed as 

oo 

S s (u) =e a ""(l + u) fcl (l-n) fe2 ^<(l + n) n , (3.1) 

n=0 

where k\ = \m — s\/2, /c 2 = |m + s\/2, and the superscript 9 denotes the association with the 
angular equation. The series coefficients are related by a three term recurrence relation and 
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the boundary condition at u = +1 is satisfied only by its minimal solution sequence. The 
three-term recurrence relation can be written as 

a e a{ + p e a 9 = 0, (3.2) 

«X+i+/3X + 7X-i = 0, (n>l), (3.3) 



where 



a e n = -2(n + l)(n + 2fci + l), 

(3 6 n = n (n - 1) + 2n(h + k 2 + 1 - 2au) - 2au(2k 1 + s + 1) 



+ (h + k 2 )(h + k 2 + l)- [a 2 oo 2 + s(s + 1) + A, 



t 'in 



7^ = 2au{n + k 1 + k 2 + s). (3.4) 

We will obtain the minimal solution if the angular separation constant A\ m is a root of the 
continued fraction equation 

ft-ft-ft-R-"" [ } 

The QNMs of the Kerr-Newman black hole are defined to be the modes with purely ingoing 
waves at the event horizon and purely outgoing waves at infinity Then, the boundary 
conditions on wave function \l/ s at the horizon (r = r+) and infinity (r — > +oo) can be 
expressed as 

1 +J + ' (3.6) 

where <7± = r+ ^ r [(r 2 + a 2 )cj — ma] . 

A solution to Eq. (j2.14|) which satisfies the desired behavior at the boundary can be written 
in the form 



q s = Vr 2 + a 2 (r - r + )-^+(r - r^ 1 ^ e Mr-r-) a n ( £± J . (3.7) 



n=0 



If we take r + + r_ = 1 and 6 = r + — r_, the sequence of the expansion coefficients {a ra : n = 
1, 2, ....} is determined by a three-term recurrence relation staring with clq = 1: 

a ai + /3 a = 0, 

a n a n+ i + + 7n«n-l = 0, (71=1,2,...)- (3.8) 

The recurrence coefficient a n , n and 7 n are given in terms of n and the black hole parameters 
by 

a n = n 2 + (C + l)n + C , 
^ = _2n 2 + (d + 2)n + C 3 , 

7n = n 2 + (C 2 - 3)n + C 4 - C 2 + 2, (3.9) 
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and the intermediate constants C n are defined by 

r, . 2i \fr 2 ,+r 2 __ „\ 
Cq = 1 — s — iuj — I — ha \ uj — ma , 




C 2 = s + 3 — 3iuj — — ( + h a 2 J uj — ma , 

C 3 = uj 2 (A + 26 - Ar + r_ + 3a 2 ) - 2mauj - s - 1 + (2 + b)iuj - A, 




lm 



uj — ma 



(3.10) 



It is interesting to note that the three-term recursion relation is obtained form the wave 
equation (|2.14jl directly. The calculation is simpler than other cases. For example, the co- 
efficients of the expansion for the electromagnetic and gravitational fields in the Reissner- 
Nordstrom black hole are determined by a four-term recursion relation which should be re- 
duced to a threaten* relation using a Lssian eliminated step (see fot detai.s). Out 
results are also more concise than that for the electromagnetic and gravitational perturbations 

n 

in the Kerr-Newman black-hole spacetime obtained by Berti and Kokkotas 25J. 

The series ()3.7j) converges and the r = +00 boundary condition (j3.6j) is satisfied if, for a 
given s, a, m and Ai m , the frequency uj is a root of the continued fraction equation 



The solution to the radial equation (|3.11|) is strictly bounded up with the determination of 
the eigenvalue A [m of the angular equation (|3.5[) . which appears explicitly in the coefficients of 
the recurrence relation associated with the radial equation. We could attempt to solve for both 
uj and A im simultaneously, but this would be a time-consuming procedure. Instead, we first fix 
the value of a, I, m and uj and find the angular separation constant A\ m by looking for the zero 
of the angular continued fraction (|3.5jl . Then, we use the corresponding eigenvalue to look for 
the zeros of the radial continued fraction. The n-th quasinormal frequency is (numerically) 
the most stable root of the n-th inversion of the continued fraction relation (J3.11|) . 

IV. NUMERICAL RESULTS 

In this section we present the numerical results obtained by using the numerical procedure 
just outlined in the previous section. The results will be organized into three subsections: 
Dirac QNMs of the Reissner-Nordstrom black hole, Dirac QNMs of the Kerr black hole and 
Dirac QNMs of the Kerr-Newman black hole. 



a n -iln a n -2ln-l «07l 
Ai-1— Pn-2— PO 



Olnln+1 a n +lln+2 « TO+2 7n+3 
Pn+2— Pn+3 - 



(n = l,2...). (3.11) 
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TABLE I: Dirac quasinormal frequencies of the Reissner-Nordstrom black hole for A = 1. 
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lo (n = 


= o) 




u> in 




= 1) 




u> (n 




2) 




uj (n 




3) 


0.00 





oooyzo- 


n 
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193965i 





295644- 





633857i 





240740- 


1 


12845i 





208512- 


1 


63397i 


0.05 





oDOoyi- 
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194066i 





296432- 





634048i 





241590- 


1 


12866i 





209400- 


1 


63420i 


0.10 





oOoOlO- 
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194366i 





298836- 





634601i 





244190- 


1 


12923i 





212121- 


1 


63480i 


0.15 





o < zuyo- 


n 
U 


19485H 





302988- 





635442i 





248696- 


1 


13001i 





216843- 


1 


635551 


0.20 





377208- 





195488i 





309132- 





636419i 





255389- 


1 


13068i 





223871- 


1 


63594i 


0.25 





384241- 





196213i 





317661- 





637242i 





264718- 


1 


13064i 





233661- 


1 


63507i 


0.30 





393651- 





19689H 





329196- 





637331i 





277339- 


1 


128711 





246805- 


1 


631181 


0.35 





406184- 





197220i 





344692- 





635457i 





294100- 


1 


12248i 





263748- 


1 


62065i 


0.40 





423108- 





196443i 





365519- 





628639i 





315375- 


1 


1061H 





282749- 


1 


59477i 


0.45 





446635- 





192202i 





391957- 





607831i 





333021- 


1 


06224i 





277400- 


1 


52872i 



A. Dirac QNMs of the Reissner-Nordstrom black hole 

As the angular momentum per unit mass tends to zero (a — > 0) the Kerr-Newman metric 



Reissner-Nordstrom metric. In this limit the angular separation constant A 



lm 



reduces to the 

A 2 = (l + |) 2 3(| ( / is the quantum number characterizing the angular distribution). Thus, 
the QNMs in this static spacetime can be obtained easily. The Dirac quasinormal frequencies 
of the Reissner-Nordstrom black hole for n = 0, 1, 2, 3 and A = 1 are listed in the table / and 
that for n = 0, 1, 2, 3 and A = 2 are listed in the table II. From the tables we find that both 
the real and imaginary parts of the quasinormal frequencies decrease as the overtone number 
increases for the fixed charge and angular quantum number (i.e., Q = const, and I = const.). 
We also find that both the real and imaginary parts of the quasinormal frequencies increases 
as the angular quantum number increases (except for the cases Q = 0.45 and n = 0) for the 
fixed charge and overtone number. 



B. Dirac QNMs of the Kerr black hole 

The Dirac quasinormal frequencies of the Kerr black hole for n — 0, 1, 2, 3 and A = 1 are 
listed in the table III and that for n — 0, 1, 2, 3 and A = 2 are listed in the table IV. We find 
that both the real and imaginary parts of the quasinormal frequencies decrease as the overtone 
number increases for the fixed angular momentum per unit mass and angular quantum number. 
We also find that both the real and imaginary parts of the quasinormal frequencies increases 
as the angular quantum number increases (except for the cases a = 0.4, 0.45 and n — 0) for 



S 



TABLE II: Dirac quasinormal frequencies of the Reissner-Nordstrom black hole for A = 2. 
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195973i 
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52082i 


0.35 





838381- 





196464i 
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04363i 





670048-1 


51544i 


0.40 





871259- 





195997i 





832478- 





600485i 





772116- 


1 


03491i 





711672-1 


49828i 


0.45 





918762- 





192368i 





882263- 





587382i 





824382- 


1 


006 17i 





760227-1 


44942i 



the fixed angular momentum per unit mass and overtone number. 

By comparing the results of the Reissner-Nordstrom and Kerr black holes, it is interesting 
to note that the real part of the quasinormal frequencies of the Reissner-Nordstrom black 
hole is greater than that of the Kerr black hole, but the imaginary part of the quasinormal 
frequencies of the Reissner-Nordstrom black hole is less than that of the Kerr black hole if 
we take the same values of the charge and the angular momentum per unit mass. That is to 
say, the decay of the Dirac fields in the charged black-hole spacetime is faster than that in 
the rotating black-hole spacetime if the values of the charge and angular momentum per unit 
mass are the same. 

C. Dirac QNMs of the Kerr-Newman black hole 

The Dirac quasinormal frequencies of the Kerr-Newman black hole for n — 0, 1, 2 and A = 1 
are shown by figure 1 and those for n — 5, 6, 7 and A = 2 are shown by figure 2. The left 
columns in the Figs. ^ and 121 describe the behavior of quasinormal frequencies in the complex 
uj plane which show that the frequencies generally move counterclockwise as the charge (the 
thick curves) or the angular momentum per unit mass (the thin curves) increases. They get a 
spiral-like shape, moving out of the Schwarzschild (Q = and a = 0) or Reissner-Nordstrom 
(0 < Q < 0.5 and a = 0) values and "looping in" towards some limiting frequencies as the 
charge and angular momentum per unit mass tend to the extremal value y/ 'a 2 + Q 2 = 1/2. For 
a given A, we observe that the number of spirals increases as the overtone number increases. 
However, for a given overtone number n, the increasing A has the effect of "unwinding" the 
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TABLE III: Dirac quasinormal frequencies of the Kerr black hole for A = 1. 
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TABLE IV: Dirac quasinormal frequencies of the Kerr black hole for A = 2. 
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spirals, as we see in the two figures that the spiral begins at n = 2 for A = 1 but starts at 
n = 6 for A = 2. 

The second and last columns in the Figs. ^ and El illustrate that the real and imaginary 
parts of the quasinormal frequencies are the functions of the angular momentum per unit mass 
a. We know from these figures that both the real and imaginary parts of the frequencies are 
oscillatory functions of the angular momentum per unit mass. The oscillation starts earlier 
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FIG. 1: The left three panels show trajectories in the complex oj plane of the first three Dirac QNMs 
of the Kerr-Newman black hole for I = 1/2 and m = 0. In each panel, the thick curve corresponds to 
the Reissner-Nordstrom black hole (a=0) which is drawn by increasing Q from zero to the extremal 
limit and the three thin curves are obtained by increasing a from zero to the extremal limit for fixed 
values of the charge Q = 0, 0.2, 0.4 (left to right). These curves show that the frequencies generally 
move counterclockwise as a (or Q) increases and the number of spirals increases as the overtone 
number increases. The other panels draw Re{oj) and Im(uj) of the quasinormal frequencies versus a 
in which the curves from bottom to top refer to Q = 0, 0.2, 0.4. These panels tell us that, for n > 2, 
both the real and imaginary parts are oscillatory functions of a and the oscillations become faster as 
the overtone number increases. 

and earlier as the overtone number n grows for a fixed A, but it begins later and later as the 
angular quantum number a increases for a fixed n. We also learn that the oscillation starts 
earlier and earlier as the charge Q increases for fixed n and A. Meanwhile, the oscillation 
becomes faster as the overtone number increases for a given A, but it becomes slower as A 
increases for a given overtone number n. 
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FIG. 2: The left three panels describe the behavior of the Dirac quasinormal frequencies of the 
Kerr-Newman black hole for n = 5, 6, 7,1 = 3/2 and m = in the complex uj plane. In each 
panel, the thick curve is drawn by increasing Q from zero to the extremal limit for a = and the 
three thin curves are obtained by increasing a from zero to the extremal limit for fixed values of the 
charge Q = 0, 0.2, 0.4 (left to right). The other panels draw Re(to) and Im{u) of the quasinormal 
frequencies versus a in which the curves from bottom to top refer to Q = 0, 0.2, 0.4. These panels tell 
us that, for n > 5, both the real and imaginary parts are oscillatory functions of a and the oscillations 
begin earlier as the charge increases. 

V. SUMMARY 

The wave equations for the Dirac fields in the Kerr-Newman black-hole spacetime are ob- 
tained by means of the Newman-Penrose formalism. The expansion coefficients of the wave 
equation which satisfy appropriate boundary conditions are directly determined by a three- 
term recurrence relation. Then, the Dirac quasinormal frequencies of the Kerr-Newman black 
hole are evaluated using continued fraction approach and the results are presented by tables 
and figures. By comparing the results of the Reissner-Nordstrom and Kerr black holes, it is 
interesting to note that the real part of the quasinormal frequencies of the Reissner-Nordstrom 
black hole is greater than that of the Kerr black hole, but the imaginary part is less than that 
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of the Kerr black hole if we take the same values of the charge and the angular momentum 
per unit mass. That is to say, the decay of the Dirac fields in the charged black-hole space- 
time is faster than that in the rotating black-hole spacetime if the values of the charge and 
angular momentum per unit mass are the same. From figures we find that the frequencies in 
the complex u plane generally move counterclockwise as the charge or the angular momentum 
per unit mass increases. They get a spiral-like shape, moving out of their Schwarzschild or 
Reissner-Nordstrom values and "looping in" towards some limiting frequencies as the charge 
and angular momentum per unit mass tend to their extremal values. For a given angular 
quantum number A, we observe that the number of spirals increases as the overtone number 
increases. However, for a given overtone number, the increasing A has the effect of "unwind- 
ing" the spirals. We also find that the real and imaginary parts of the quasinormal frequencies 
are oscillatory functions of the angular momentum per unit mass, and the oscillating behav- 
ior starts earlier and earlier as the overtone number grows for a fixed A (or as the charge Q 
increases for fixed n and A), but it begins later and later as A increases for a fixed overtone 
number. 
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